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Abstract. Bivariate thematic maps are powerful tools, making visible spa- 
tial associations between two geographic phenomena. But bivariate themat- 
ic maps are inherently more visually complex than a univariate map, a 
source of frustrati on for both map creators and map readers. Despite a vari- 
ety of visual solutions for bivariate mapping, there exists few empirically 
based 'best practi ces' for sel ecti ng or i mpl ementi ng an appropri ate bi vari ate 
map type for a given scenario. 

This research reports on a controlled experiment informed by the theory of 
selective attention, a concept describing the human capacity to ignore un- 
wanted stimuli, and focus specifically to the information desired. The goal 
of this research was to examine if and how the perceptual characteristics of 
bivariate map types impact the ability of map readers to extract information 
from different bi vari ate map types. 

55 participants completed a controlled experiment in which they had to an- 
swer close ended questions using bivariate maps. The results of this exper- 
iment suggest that 1} despite longstanding hesitations regarding the utility 
of bivariate maps, participants were largely successful in extracting infor- 
mation from most if not all of the tested map types, and 2) the eight map 
types expressed unique differences in their support for different map read- 
ing tasks. Though selective attention theory could explain some of these dif- 
ferences, the performance of the map types differed appreciably from si mi- 
lar studies that examined these map symbols in a more abstract, percep- 
tion-focused setting. While the perceptual models of selective attention can 
still be useful in guiding map design, more work is needed to understand 
the cognitive aspects and limitations to bivariate map reading. 

Keywords: Bivariate Maps, Map Perception, Thematic Mapping, Selective 
Attention, User Studies 



Introduction 



Existing Thematic Cartography literature primarily focuses on univariate 
maps, or cartographic representations that portray only one attribute of a 
geographic information set. Displaying two or more attributes (a bivariate 
and multivariate map, respectively) is a powerful way to convey information 
about associated geographic phenomena, but successfully designing bivari- 
ate/ multivariate maps is challenging due to the added density of infor- 
mation. 

The functional purpose of a bivariate map is to show relationships among 
two geographic phenomena (Fisher, 1982; Tyner, 2010). Visualizing this 
geographic relationship frequently provides insight into understanding the 
mapped phenomena. The cost of bivariate mapping is in comprehensibility. 
A bivariate map is, by its nature, more visually complex than a univariate 
map. Visual complexity makes the map more difficult for the viewer to pro- 
cess mentally, and, if this complexity proves overwhelming, it can render 
the map valueless to the reader. Fisher (1982: 268) puts the multivariate 
map in simple cost/ benefit terms: there is "a limit beyond which the diffi- 
culty of comprehending two or more subjects exceeds the value of being 
able to relate them". 

This writing supposes, along with Carswell & Wickens (1990) and Nelson 
(1999; 2000) that the challenges of designing and reading bivariate maps 
can be resolved, at least i n part, by a better understandi ng of the perceptual 
and functional aspects of various bivariate mapping solutions. The parcity 
of user-based, comparative testing of bivariate map types presents a signifi- 
cant problem for the thematic map designer; with little theoretical or em- 
pirical basis from which to draw, the designer cannot make an informed 
decision on which bivariate symbolization scheme to use. There are many 
more design possibilities for bivariate mapping compared to univariate 
mapping; numerous bivariate map types are already in standard use, and 
new solutions continue to be proposed. Given this, the need to understand 
the capabilities and limitations of these various bivariate map types be- 
comes al I the more pressi ng. 



Literature Review 

2.1 Taxonomy of Bivariate Map Types 

Cartography textbooks discuss bivariate and multivariate thematic mapping 
less systematically than, for example, univariate map types, map projec- 
tions, or color schemes. An examination of different thematic cartography 
textbooks (Dent, Torguson, & Hodler 2009; Fisher, 1982; Krygier & Wood, 



2011; Slocum et al., 2003; Tyner, 2010) reveals that these sources identify 
eleven different bivariate map types, which represent only a small selection 
of map types identifiable from published maps or the cartographic litera- 
ture. No text describes more than four different bivariate map types, and 
only two map types (bivariate choropleths and choropleth with overlaid 
graduated symbols) are considered by more than two sources. 

I n the i nterest of creati ng a more exhausti ve taxonomy of bi vari ate map 
types, th i s wr i ti ng consi ders al I bi vari ate map sol uti ons to be a combi nati on 
of two visual variables (size, shape, value, etc, as initially described by Ber- 
tin, 1967|1983) and two symbol dimensionalities (point, line, polygon, etc) 
(Figure 1). This issimilar to the combinatorial systems proposed by Stefan 
et al., 2007, and Nelson, 2000. 
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Figure 1 The constituents of abivariate map symbol: two visual variables 
combined with two symbol dimensionalities. 

Different authors have offered different collections of visual variables 
(Tyner, 2010), and these listings are sure to adapt as new visual variables 
are proposed. Depending on the number of visual variables considered, the 
number of two- vari able combi nati ons can range in the dozens to hundreds. 
Symbol dimensionalities, comparatively, are generally agreed to be limited 
to points (Od), lines (Id), and polygons (2d), with other authors also consid- 
ering surfaces (2.5d) and volumes (3d) (Tyner, 2010, Stefan et al., 2007). 
The solution space for bivariate map symbols, then, is exponentially large. 
Fortunately, the number of tenable graphical solutions for a given use case 
can be minimized by considering the appropriateness of a given visual vari- 
able or symbol dimensionality to represent the data set (per MacEachren, 
1995). As an example, Figure 2 provides an example of graphic solutions 
appropriate to representing two sets of quantitative information enumerat- 
ed to polygonal features. 





Shaded Cartogram 



BivaTiate 
Choropleth 




VBA Cartogram VBA Choropleth 




Cartogram w/ 
Texture 




Cartogram w/ 
Dot Density 

Size 




Shaded Texture 



Shaded Dot 
Density 



VBA w/ Texture 



VBA Dot 
Density 



Polygon/Polygon 
Solutions: 

Both visual variables 
are applied to the 
map features. 



Empi/SBmsiBpfaseti sofcdbns tataro 
rxfrhMXtoxietorgraphkegyunianatih. 
Map typos in .Safes demtos aoUnns 
Ink 




Bivariate Tejrture 



Color Transparency 





FUlDrmstly 




Cartogram w/ 
Graduated Symbols 




Choropleth w/ 
Graduated Symbols 




o 



VBA w/ Graduated 
Symbols 




Cartogram w/ 
Rotated Symbols 




Choropleth wf 
Shaded Symbols 




O 



VBA w/ Shade] 
Symbols 



Choropleth wf 
Rotated Spobdls 



VBA w/ VBA 
Symbols 



P51 




VBAw/ 
Rotated Symbols 



Point/Polygon 
Solutions: 

One visual variable is applied 
to the map features, one is 
applied to a point symbol. 



Empty tamaiBprBaenteaiuSkins that are 





Texture w/ 
Graduated Symbols 



Texture w/ 
Shaded Symbols 




Texture w/ 
VBA Symbols 




Texture w/ 
Rotated Symbols 



Dot Density w/ 
Graduated Symbols 





Dnt Density w/ 
Shaded Symbols. 



Dot Density w/ 
VBA Symbols 




Dnt Density w/ 
Rotated Symbols 



Color 



Tnanapaaxneg Orientation 



FDlDrawiliy 



i 




IB 




tfectana/e Map 



Split Proportional 
Symbol (size) 



To] 



o 




Graduated Shaded Split Proportional 
Symbols Symbol (color) 



Graduated VBA 
Symbols 





Graduated Shaded Rotated 

Rotated Symbols Symbols 



Point/Point 
Solutions: 

Both visual variables are 
applied to a point symbol. 



Empty iteim reprasert saltations that aw 




Split Proportional 
Symbol CVBA) 



"0T"! 



Rotated VBA 
Symbols 




Color Transparency Orientation FBlSfee niLBcakrfty 

Figure 2. A taxonomy of bivari ate symbol solutions appropriate for two sets of 
numerical data aggregated to a polygonal enumeration unit. Hue, saturation, and 
value have been combined into simply 'color' to avoid redundancy. 



2.2 Bivariate Maps and Selective Attention 

Defined simply, selective attention is the ability of an observer to attend to 
one stimulus while ignoring the confounding influence of others. In the 
context of complex visual stimuli (such as a map), selective attention 
manifests as the ability to attend to specific visual variables while ignoring 
the others. The theory of selective attention was developed by research 
psychologists, and has been proposed as a means to aid understanding and 
design of information graphics, including thematic maps (Carswell & 
Wickens, 1990; Shortridge, 1982; Nelson, 1999;2000). One contribution 
from Selective Attention research has been the concepts of separability, 
integrality, configurality, and asymmetry (ibid). These phenomena arise 
when two visual variables are combined in a symbol, describing the 
perceptual aspects of that visual combination. 

In a separable combination, the viewer is capable of attending each 
individual visual variable, with minimal perceptual i nterference f rom the 
other. I n an integral combination, the visual variables interact in such a way 
that attention to the individual variables are inhibited. In Figure 3, for 
instance, notice how the 'cluster' in the lower left is easier to notice as four 
squares in the separable combination, as opposed to four rectangles of 
equal width. Contrarily, it is easier to spot the rectangles with a large area 
as opposed to dark hexagons and pentagons. 

Configurality exhibits properties of both separability and integrality, and 
tend to arise when the symbol employs the same visual variable twice. 
Uniquely, configural combinations have features that arise when the 
variables are 'in agreement', such as the unintereupted circles in Figure 4. 
Finally, asymmetrical combinations appear to represent unique or 
conditional interaction effects between the two visual variables. I n Figure 4, 
for instance, variations in hue are least noti cable when the symbols are very 
dark. 
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Figure 3. The same data represented using a separable combination (shape/ value, 
left) and an integral combination (height/ width, right). 
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Figure 4. A configural combination (value/ value, left) and an asymmetrical 



combination (hue/ value, right). 

2.3 Bivariate Maps & Encoding Information 

I n a bivariate map symbol, the two statistical attributes are encoded by the 
two visual variables used to construct the symbol. Each combination of vis- 
ual variables also creates two emergent (gestalt) visual dimensions, which 
are visual variables unto themselves. A combination of height/width, for 
instance, creates two emergent visual dimensions: the total area of the 
symbol, and its directionality (vertical vs horizontal). 

These gestalt visual variables encode information in the same way the orig- 
inal visual variables do. One encodes a positive association between the da- 
ta (i.e., map features where both attributes are low vs features where both 
attributes are high) whilethe other gestalt dimension encodes negative as- 
sociation (features where one data variable is higher than the other, or vice 
versa). Given a matrix of symbols (as would be seen in a two-axis map leg- 
end), these emergent visual variables can be found along the orthogonal 
axes (Figure 5). (These orthogonal data axes will be referred to as the 
Plus(+) and Minus(-) axes, and the original information axesastheX and Y 
axes). 
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Figure 5. The orthogonal data axes (above), and two examples of the gestalt visual 
variables created within a bivariate map symbol. 



Methods 

A controlled experiment was administered to assess the variation in per- 
formance across various bivariate map design solutions. The purpose of the 
study was to examine empirically how different map types and different 
conditions of selectivity support different map readi ng tasks, and if this lev- 
el of support varies according to differences in the expertise of the map 
reader. Conceptually, the experiment was designed to reconcile, at least in 
part, the methods employed in Selective Attention research and the meth- 
ods employed in cartographic performance testing. This involved asking 
parti ci pants to perform a selection of map readi ng tasks across eight differ- 
ent bi vari ate map types, recordi ng thei r accuracy and response ti me to each 
question. The survey opened with questions designed to assess each partic- 
ipants' expertise with regards to map use, and closed with questions on the 



users' personal preferences of the different map types, modified from ques- 
tions used by Olson (1982). 



A total of 55 participants participated in the study. The majority of partici- 
pants were recruited from the Geography Department of the University of 
Wisconsin - Madison, although the study was open to any i nterested partic- 
ipants. Subjects were recruited purposefully to represent a range of experi- 
ence and knowledge of cartography and spatial analysis. The study partici- 
pants were offered 5 U SD as compensati on for thei r ti me. 
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Figure 6. Legends for the eight map types tested. 

Eight common bivariate map types, drawn from the map types discussed in 
Figure 2, were designed for inclusion in the experiment (Figure 6). Two 
map types were included for each condition of selectivity (separable, inte- 
gral, configural, and asymmetrical). The map types were selected based on 
the foil owing criteria: 1) if possible, there was existing research establishing 
the visual combination as separable, integral, etc., 2) the map type was par- 
ticularly representative of the perceptual features associated with that selec- 
tivity, and 3) the map type is commonly used in cartographic practice. 

Note that the design for the bivariate choropleth uses a modified construc- 
tion from the other legend designs: the visual variables are applied across 
the orthogonal information axes, and it employs small variations in satura- 
tion as well as hue. This was done to align the legend design with estab- 
lished recommendations for the design of bivariate choropleth legends 
(Trumbo, 1981; Dunn 1989). The design for the value-by-alpha legend in- 
cluded a pattern fill to allow for a better depiction of transparency; on a 



matte white or black background, variations in transparency are indistin- 
guishablefrom variations in value, rendering a value- by- alpha map concep- 
tually identical to a bivariate choropleth map (Roth, Woodruff, &J ohnson 
2010). 

Each map legend was applied to a basemap composed of 36 counties from 
western Ohio, rotated 90 degrees, and modified slightly in size and topolo- 
gy. The allotment of data values within each map was normalized and ran- 
domized. Each map included a legend depicting the fictitious variables of 
'pizza consumption' and 'chicken consumption' along an ordinal scale 
(high/ medium/low) (Figure 7). 
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Figure 7. Screenshot from the test interface. 

Each map type had eight accompanying questions, requiring participants to 
compare the values of map features across two levels of reading across and 
four search axes. The levels of reading included Elementary (looking at sin- 
gle map units) and General (looking at regions composed of multiple units). 
The search axes included the information axes described in Section 2.3: X 
(the first mapped attribute), Y (the second mapped attribute), Plus(+) (pos- 
itive association between attributes; High/High to Low/Low), and Minus(-) 
(negative association between attributes; High/Low to Low/High). Each of 



these 'blocks' of eight questions opened with a screen that permitted partic- 
ipants to examine the example map and legend. The order of map types 
shown to participants, as well as the order of questions within each block, 
was randomized. 



Results 

3.1 Statistical Results 

Participants were consistently successful in accurately answering the 64 
questions (Table 1). The participants' global accuracy rate was 96.1%, and 
26% of the map questions were answered correctly by every single partici- 
pant. ANOVA analysis of accuracy rates between trials did not show statisti- 
cally significant differences between accuracy across the different map 
types (p- value = 0.064) nor across the different conditions of selectivity (p- 
value = 0.0598). 

Response times (RTs) showed statistically significant differences across 
both map type, map selectivity, and task (P <2e-16, using ANOVA) (Figure 
8). These differences indicate that the choice of map type has important 
impacts on how intuitively users are ableto extract various forms of infor- 
mation from the map. The best performing maps (that is, those with the 
lowest average RT combining all tasks) were the rectangle map (2115 sees) 
and the choropleth with graduated symbols (22.97 sees). The worst per- 
formers were the value- by-alpha (29.39 sees) and spokeglyph (38.33 sees). 
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Table L Accuracy rates by map type and task. Table divided for legibility. 




Figure 8. Reaction Times (RTs) by map type and task. Left side of bar represents 
elementary- level task, right side represents general level. 



The closing questions on user preference showed clear differences among 
the eight map types (Table 2). Kruskal-Wallis testing (a non- parametric 
version of ANOVA) found statistically significant differences between how 
users ranked the map types (p value = 2.2e-16). Table 3 provides a sum- 
mary of the Likert scores. There is little variation in how participants re- 
sponded to the questions within each map type, however: that is, if a partic- 
ipant liked or disliked a map type, they would rate it similarly high or low 
regardl ess of the questi on asked. 

Participants' accuracy and response times did not significantly vary accord- 
ing to user expertise (Table 3). Participants reporting an educational back- 
ground in cartography/ GIS had an overall accuracy merely 0.04% better 
than those without. Nearly identical results occur using the other major 
measure of expertise, work experience. Again, no statistically significant 
differences were found between those with and those without work experi- 
ence i n Cartography/ Gl S. 
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Table 2. Likert scale results: mean and standard deviation for participants' scoring 
of the map types along various scales (all running 1-7, with 1 representing the de- 
scriptor on the left, and 7 the descriptor on the right). 
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Table 3. Global Response time (in seconds) and accuracy rates across cartog- 
raphy/GI S related user expertise. 

3.2 Analysis 

Despite hesitations about the utility of bivariate maps (described in Section 
1), the participants in this survey were consistently successful in answering 
the questions presented to them. There were also several map types (such 
as the choropleth w/ graduated symbols and bivariate choropleth) that us- 
ers rated largely positively in their capacity to read and understand the in- 
formation on the map (based on the Likert scores). It should be recognized, 
however, that the structured, task- specific way participants interacted with 
the maps in this study would have engaged different visual and cognitive 
activities than the unstructured, spontaneous kinds of map reading that is 
more I i kely to occur i n everyday I ife (Antes & Chang 1990). 

Selective attention theory provided some insights into the performance of 
the eight map types, but could not account for all results seen. The best per- 
formance in the General Plus(+) task, for instance, was expected from the 
strong emergent dimensions of the integral map types, however the best 
performer in this task was the bar chart, a configural solution. The X and Y 
tasks were hypothesized to be best supported by the separable combina- 
tions, but this materialized in only one instance (the choropleth with gradu- 
ated symbols, in the General Y task). Additionally, participants did not rate 
either separable or integral combinations as better at portraying individual 
distributions versus relational information. A given map type did not al- 
ways perform identically to the other map type with the same selectivity; 
the disparities in performance between the configural solutions (bar chart 
and spoke glyph) provide a good example. 

The eight map types showed frequent variations in reaction times across 
task, whether that be performing better than the other maps in one of the 
eight tasks, or particular map types showing variation in performance be- 
tween Elementary and General level tasks within a given information axis 



(Figure 8). There were also statistically significant differences between the 
eight tasks in terms of their accuracy rates. Speaking broadly, variations in 
task performance across the map types were rel atabl e to the uni que percep- 
tual properties of their specific visual cues. In one case, the spoke glyph, the 
lack of strong visual cues made it challenging for participants to retrieve 
information regardless of task. The other map types were generally success- 
ful in supporting the eight tasks overall, but showed differences in perfor- 
mance dependent on the information axis or the level of reading. The con- 
figural map types, for instance, demonstrated a good ability for viewers to 
extract information along the Mi nus(-) axis, but only when attending to in- 
dividual symbols (Elementary tasks); when challenged to perform the same 
task at the General level, participants seemed unable to visually aggregate 
the symbol s i nto a si mi I ar I y i ntui ti ve cue. 



Conclusions 

The major goals of this research were to 1) understand the unique percep- 
tual qualities of different bivariate mapping solutions and 2) empirically 
test if these qualities impacted various forms of information retrieval. Alt- 
hough selective attention provided a useful framework to understand the 
perceptual characteristics of the map types, the results of the study differed 
from expectations derived from earlier research, which examined similar 
symbols in an abstracted setting. 

This study assumed that the major challenges of reading a bivariate map 
were perceptual : that is, the visual variables interact with each other in such 
a way as to hinder information extraction from the map, and choice of map 
type could serve to minimize visual interference and maximize the intuitive 
representation of the information (or, at least, certain aspects of the infor- 
mation). But cognition appears to also be an important aspect of bivariate 
map reading. Other than the Likert scales, the information collected in the 
survey demonstrates very little about how the maps were cognitively under- 
stood by the participants. I n order to understand why reaction times varied 
across task, for instance, it would have been beneficial to know the mental 
strategies the participants used to answer the questions. Did the partici- 
pants begi n by looki ng at the map, or at the question itself? H ow often, and 
in what contexts, did participants have to re-check the legend? Did they an- 
swer the General -level questions by attending to individual map features, or 
by trying to visually aggregate them? There are several methods that could 
better investigate such cognitive- based questions, such as eye tracking, fo- 
cus groups, "think-aloud" experiments (Pickle 2003), or by giving parti ci- 



pants more open-ended information-seeking challenges (rather than the 
highly specific information retrieval tasks used in thisstudy). 

Similarly, this research attempted to identify and discuss the emergent di- 
mensions and features of various bivariate map types; dimensions such as 
the 'leftedness/rightedness' of a split proportional symbol, or the direction- 
ality (vertical/ horizontal) of a rectangle map. The understanding of these 
various emergent visual dimensions is significantly less than the under- 
standing of the more fundamental visual variables. Since these emergent 
dimensions do serve to encode information in a bivariate map, it would be 
enlightening to examine how these dimensions are perceived and their ap- 
propriateness for encoding different varieties of data. 
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